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SUMMARY 


Interferograms were taken of the flow around a sphere at Mach num- 
bers of 1.30 and 1 . 62 . The results of the evaluation of these inter- 
ferograms are shown as plots of contours of constant density ratio in 
the flow field around the sphere and as the distribution of flow vari- 
ables along the axis of symmetry between the shock wave and the sphere . 
In appendixes the theory of the analysis is reviewed, sources of error 
are discussed, and the detailed procedure used in the evaluation of the 
interferograms is described. 


INTRODUCTION 


One of the flow problems currently receiving attention is that of 
the flow behind detached shock waves. The present paper gives experi- 
mental data on the flow behind the detached shock waves on a sphere. 
These data consist principally of density distributions around spheres 
and were obtained for two supersonic Mach numbers by analysis of inter- 
ferograms. Previous investigations of axially symmetric flow by inter- 
ferometry can be found in reference 1. 

For the convenience of those interested in the application of 
interferometry to the study of axially symmetric flow fields, the 
symbols used are given in appendix A and the theory of the analysis is 
reviewed in appendix B. A discussion is given in appendix C of the 
sources of error and the accuracy of the results obtainable from inter- 
ferograms of axially symmetric flow. (This discussion is somewhat 
similar to those given in references 1 and 2.) Because reasonable 
accuracy in the results requires that much care be taken in measuring 
and in plotting fringe shifts, a detailed description of the computa- 
tional procedure that was used to obtain the final results is given in 
appendix D. A table of the coefficients for 50 zones for use in the 
Weyl method of evaluation has been calculated and is also given. 
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APPARATUS 


The test setup consisted principally of a supersonic- jet apparatus 
that gave a 3-inch-square free jet and an interferometer with 4-inch- 
square plates. This apparatus has been described in reference 3- The 
only change in the apparatus was a change in the film magazine to permit 
the use of 70 -fliillimeter film rather than the 35-millimeter size previ- 
ously used. Two different Mach numbers were obtained by using two 
supersonic nozzles similar to those described in reference 3- The test 
model was a steel sphere l/2 inch in diameter that was supported from 
behind by a sting l/8 inch in diameter. 


RESULTS AND DISCUSSION 
Interferograms 


An interferogram that showed the undisturbed fringes in the test 
section was taken with no flow (fig. 1) and then an interferogram of 
the flow around the sphere at Mach number 1.30 was taken (fig. 2 ). The 
same procedure was followed at a flow Mach number of 1 . 62 ; the no-flow 
interferogram is shown in figure 3 and the flow interferogram, in fig- 
ure 4. Cross sections were drawn on the flow interferograms perpendicu- 
lar to the axis of symmetry. The variation of density along each cross 
section was then obtained by the method described in appendixes B and D. 
A plot of the ratio of density along the cross section to free-stream 
density was made for each cross section. 


Density Distributions 

The plots of density ratio along the cross sections were used for 
plotting contours of constant density ratio in the flow field around 
the sphere. Experimental points taken from the unfaired plots of den- 
sity ratio are shown in figure 5* The curves that are shown are faired 
through these points. 

The method of analysis of the interferograms that was used has the 
disadvantage that the results obtained are inaccurate in a small region 
immediately behind a discontinuity in the density. (The reason for the 
inaccuracy is explained in appendix C.) Because of this inaccuracy in 
the results imm ediately behind the shock wave, points on the plots of 
density ratio that were near the shock wave and that were obviously in 
error were not included in figure 5. For each contour in figure 5> 
therefore, a small region adjacent to the shock wave exists in which 
experimental points have not been obtained. The density contours have 
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not been extended through this region to the shock wave. On the shock 
wave , however , the correct location of the intersection of each contour 
with the shock wave has been indicated. The locations of the inter- 
sections were obtained from the equation for the density ratio across a 
two-dimensional oblique shock wave with use of measured slopes of the 
shock wave. 

The magnitude of any systematic errors in the contours shown in 
figure 5 is not known. The uncertainty in the location of the contours, 
aside from any systematic errors, is believed to be not larger than plus 
or minus approximately one-half the distance between adjacent contours. 


Flow Variables along Axis 

The distribution of the flow variables along the axis of symmetry 
is of interest. The region of the axis is, however, a region of rela- 
tively large error. In order to obtain what are thought to be the best 
values of density ratio along the axis, the values at the axis, obtained 
by fairing the plots of density ratio for each cross section, were used 
to plot the variation of density ratio along the axis between the shock 
wave and the sphere (fig. 6 ). Immediately behind the normal part of the 
shock wave, the density ratio is well-known from the Rankine-Hugoniot 
relations. At the nose of the sphere, the correct value of the stagna- 
tion density is also known. These theoretical values were used in 
fairing the curves shown in figure 6 through the experimental points. 

The density ratios along the axis given by the faired curves have also 
been converted to ratio of static pressure to stagnation pressure 
(fig- 7). The distribution of density ratio along the axis has also 
been converted into the distribution of Mach number shown in figure 8 . 
(The stagnation temperature for the free-stream Mach number of 1.30 
was 524. k° F abs. and for Mach number 1.62 was 533.1° F abs. The free- 
stream density for Mach number 1.30 was 0.003215 slug/cu ft and for 
Mach number 1.62 was 0.003501 slug/cu ft.) 

The variation of density ratio along the surface of the sphere for 
Mach number 1.62 is shown in figure 9 . The experimental points were 
taken from the plots of density ratio along each cross section. 


CONCLUDING REMARKS 


Experimental data on the flow around spheres at supersonic speeds 
have been presented. Although greater accuracy in the data would be 
desirable, the data given herein are accurate enough that they should 
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be of use in checking and evaluating theoretical treatments of the 
problem of predicting the flow behind detached shock waves. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., June 12, 1950 
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APPENDIX A 
SYMBOLS 


a 

A k,i 

d 

i 

k 

L 

M 

N 

n 

P 

r 


R 

S 

T 



w 

X 

y 


z 

7 


radius of first zone 
coefficient 

axial distance from center of curvature of shock wave 
number of zone at which density is required 
zone number 

distance through test section or through disturbance 
Mach number 

number of outermost zone 
index of refraction 
pressure 

radial distance from axis of symmetry 
radius 

fringe shift per unit fringe width 
temperature 


value of v for which density is required 

Cartesian coordinate in direction of undeviated light ray 

Cartesian coordinate perpendicular to light ray and to free- 
stream flow direction 

zone width 

ratio of specific heats 
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K Gladstone -Dale constant 

X 0 wavelength of light in vacuum 

v difference between index of refraction at point r 

reference index of refraction 

p density 

p reference density 

ref 

0 angle between ray and normal to surface 

Subscripts: 
s stagnation 

0 free stream 

1 unrefracted rayj also, correct value 

2 refracted ray; also, approximation to correct value 


and 
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APPENDIX B 
ANALYSIS 


The problem is to determine the density from the known fringe 
shift. For the two-dimensional case, where the density is not a func- 
tion of x, the relation between fringe shift and density has been 
shown in reference 3, for example, to be 



For the three-dimensional case, where the density is a function also 
of x, the fringe shift at any point on the interferogram is the sum- 
mation of the shifts produced in each element of path of the ray of 
light. If the assumption is made that the density distribution is 
axially symmetric, then p is a function of r alone. If the addi- 
tional assumption is made that, within the disturbance, the light ray 
does not deviate from a straight line that is parallel to the x-axis, 
then the relation between fringe shift and density is 



Three methods of handling equation (l) have been proposed. The 
earliest is that which is discussed by Schardin (reference 4 ). In this 
method, the disturbance is divided into zones, as shown in figure 10, 
where a zone is the region between two adjacent concentric circles. 

(The zone between circles i and i + 1, for example, is referred to 
in the present analysis as zone i.) The assumption is then made that 
the density is constant in each zone. The resulting equation, which is 
thus an approximation to equation (l), can be solved for p(r). 

In the second method, which was introduced by Van Voorhis (refer- 
ence l), the value of the quantity p(r) - P re p in each zone is assumed 

to be a linear function of the radius r of the zone. The resulting 
equation, which is thus also an approximation to equation (l), is 
solvable . 


In the third method, which was introduced by Weyl (reference 2), 
equation (l) is solved exactly for p(r) - p ref and the resulting 

definite integral is approximated by assuming that in each zone the 
fringe shift is a linear function of y^. First the variables 
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v = 


2 


r 


u = 


and 



are substituted into equation (l) to give the equation 

s(u> - i p-MJ? 

K oJ u \jv — u 

This integral equation is the same as Volterra's integral equation of 
the first kind except that the limits of integration differ from the 
Yolterra equation limits of 0 and u. This equation can be solved 
in a manner similar to that used for solving Volterra's equation (refer- 
ences 2 and 5). The result is 


v(w) 



S 1 (u) du 
\ju - w 


where the prime indicates differentiation with respect to u. This 
resulting integral can be computed by breaking the integral into steps, 
one step for each zone, and by assuming that in each zone the fringe 

shift S(u) is a linear function of u, where = u = r i+l 2- T^ en 


v i = 


r U i+l S 1 (u) du 
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For zones of equal width z y 


and 


r 


i 


iz 


V . = 


2 V& 

7 . Z- 


TtZ 


k=i 


n \ \/(k + l) 2 - i 2 - \/k 2 - i 2 

S k+l) 2 F+ 1 


( 2 ) 


Equation (2) is an approximation to the solution of- equation (l) and 
is the relation that was used for computing the results given in the 
present paper. 


The coefficients 


kj i 


. _ \/(k + i) 2 - i 2 - Vi 


k 2 - i 2 


2k + 1 


for 50 zones were calculated by the Bell Telephone Laboratories 
X-667LL relay computer at the Langley Laboratory and are given in 
table I . 


The situation may arise that would make desirable the use of more 
than 50 zones. Because the number of coefficients is approximately 
equal to .half the square of the number of zone s , a large number of coef- 
ficients would have to be calculated. A method has been devised in 
which the number of coefficients is equal only to the number of zones 
and each coefficient is also easier to calculate. The idea is to use 
zones not of equal width but of radius proportional to the square root 
of i. 


r = a \fT 

z = r i+l - r i 

= a(\/i + 1 - x/T) 

where a is the radius of the first zone. Then 


. . / , , \ 

i ' ifL ( s k - ®k+l) (v/kU 1 + 1 -\/k - ij 

k=i ' 
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The coefficient is a function only of k - i, which ranges in value 
from 0 to N - 1. Therefore, only W coefficients need be calculated, 
rather than N^/2. This new method was not used in the calculations for 
the present paper but is mentioned because of its possible value in 
future calculations. 
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APPENDIX C 
SOURCES OF ERROR 


The purpose of the present section is to consider some of the 
sources of error inherent in the application of equation (2). Equa- 
tion (2) can be used for determining density distribution with suffi- 
cient accuracy only if the following conditions are suf f iciently well 
satisfied: 

(a) The index of refraction of the gas varies linearly with density 

(b) A reference density at some point in the field of view of the 
interferometer can be determined 

(c) The disturbance is rotationally symmetrical about an axis that 
is known 

(d) The light beam is not refracted 

(e) The variation with y of fringe shift S(y) for a cross 
section of the disturbance perpendicular to the axis of symmetry is 
known 

(f) The variation with u of fringe shift S(u) can be replaced 
by a series of straight lines 

As far as is known, condition (a) is always sufficiently well 
satisfied by air. Whether condition (b) is satisfied at all is gener- 
ally determined by the experimental setup. For the present investiga- 
tion the reference density was the free-stream density. For determining 
the free-stream density, the value of the stagnation temperature was 
measured at the instant the interferogram was taken. The Mach number 
given by the nozzle was known from a previous calibration described in 
reference 6. The free-stream temperature was then calculated from the 
relation 



M 2 


Then the free-stream density was computed by the general gas law under 
the assumption that the free-stream static pressure was equal to the 
ambient atmospheric pressure, which was measured. The determination of 
free-stream density is, of course, subject to some degree of uncertainty. 
The error in free-stream density is estimated to be not greater than 
3 percent, and the probable error is estimated to be about 1 percent. 
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Condition (c) is assumed to be satisfied for the part of the flow 
field that is behind the shock wave and ahead of the most upstream part 
of the intersection of the shock wave with the mixing region at the 
edges of the jet. Because the shock wave is axially symmetric and the 
jet is square, as soon as the shock wave reaches the edge of the jet, 
condition (c) is no longer satisfied. 

Condition (d), that the light beam is not refracted, is usually 
sufficiently well satisfied in rotationally symmetric flow. (The only 
place the effects of refraction are likely to be large is in boundary 
layers in two-dimensional flow. ) That the effect of refraction on 
fringe shift and on calculated density was fairly small in the present 
case was determined by calculating the effect. The purpose of the 
calculation was not to make a precise analysis of the effects of refrac- 
tion but to determine approximately the magnitude of the effect on cal- 
culated density in the region where the effect was largest. That region 
is at the nose of the shock wave. Certain assumptions that greatly 
simplified the analysis were therefore justified. The assumption was 
made that the region behind the nose of the shock wave was a spherical 
region of constant density. When a ray of light impinges obliquely on 
a region of greater density, the ray is refracted and emerges at a 
place different from the place it would emerge if not refracted. The 
optical path length is, therefore, changed by the refraction. Because 
of the change in optical path, the fringe shift produced by the ray is 
also changed. Figure 11 shows the ray of light entering the sphere 
at Pi and leaving at P 2 . Only a meridian plane, defined by the ray 
and the center of the sphere, need be considered. The geometrical path 
length without refraction is Li and with refraction is Lq 


Lq = 2R cos 0q 


= 2R\jl - sin 2 0 x 


L 2 = 2R cos 02 



sin 0i = | 


and, inasmuch as 


sin 0 2 = ^ sin 0.^ 
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then 


sin 02 


npd 

r^R 


Therefore, 



The value of np outside the sphere (in the free stream) is 1.0004 at 
Mach number 1.62. The density within the sphere is assumed to be twice 
that outside the sphere. Then, because n - 1 is proportional to the 
density, the value of n2 inside the sphere is 1.0008. The ratio of np 
to np is 1.0004. The values of I^/Lp for various values of. d/R 
are shown in the following table, together with the percentage error in 
the fringe shift. The corresponding values of distance behind the nose 
of the assumed sphere of constant density, expressed in terms of the 
distance between the nose of the shock wave and the nose of the sphere, 
are shown in the last column of the table. 


d/R 


Error in 
fringe shift, 
percent 

Percent of 
distance from shock 
wave to sphere 

0.995 

1.039 

3.9 

3 

.990 

1.019 

1.9 

5 

.980 

1.010 

1.0 

11 

.960 

1.005 

.5 

22 

.9^0 

1.003 

.3 

32 

.920 

1.002 

.2 

^3 


The values of the radius of curvature of the shock wave at its nose, 

2.9 times the radius of the model, and of the axial distance between 
the shock wave and the model, O.56 times the radius of the model, were 
taken from reference 6. 

The effect of refraction is seen to be small except immediately 
behind the shock wave. There, however, the error due to refraction may 
still be relatively unimportant because at the same place there may be 
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a larger error due to the "overshooting" in the computed density that 
occurs immediately behind a discontinuity in the density, such as exists 
at a shock wave. (The overshooting is explained subsequently in the 
discussion of condition (f).) 

Whether condition (e) is sufficiently well satisfied for the com- 
puted value of density to lie within given limits of uncertainty depends 
on a number of factors. A relation that gives the effect of some of 
these factors is derived. Figure 10 shows a cross section of the flow 
field divided into concentric zones. A typical ray of light is shown 
traversing the field. The innermost zone through which this ray travels 
is designated zone i. The assumption is then made that, for this ray, 
the fringe shift due to the density outside of zone i is accurately 
known. The problem is to determine how much error in the fringe shift 
at y = i can be tolerated for a given uncertainty in the determination 
of the density at the same place. Let Sj^ be the correct value of 

fringe shift at y = i and Sj^ the approximation to S-^. The error 

in fringe shift is 


AS - Sj_^ - Sj^ 

From the definition of v and equation (2), 

V i x = - P T ef) 


" *z Aj- ( Sk ' Sk+1 ) 1 Ak >i 


v i 2 K v 2 " p ref 


Jtz 


N-l 


k=i 



3 k+l 


2 A k,i 


Then, because all the fringe shifts from k = i + 1 to k = N - 1 are 
assumed to be accurately known and are the same for both case 1 and 
case 2, 


Ap = p 1 - p 2 


2K 


AS A i>jL 


TUZk 


(3) 
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Equation (3) can be used to calculate the accuracy with which fringe 
shift must be known for determining density within specified limits. 

The equation is- used herein to investigate the effect of error in 
the value of the fringe shift on the value of the density for a sample 
case. A square density distribution as shown in figure 12 is assumed. 
The corresponding fringe-shift curve is also shown. The assumption is 
made that this fringe-shift curve is known to ±0.001. From this fringe- 
shift curve the density distribution is in turn recalculated. The 
result is the solid curve that peaks at zone 24. The overshooting that 
occurs at the discontinuity in density is not the result of the uncer- 
tainty of 0.001 in the fringe shift at each zone but is caused by fail- 
ure of condition (f) to be satisfied and is discussed subsequently. 

For comparison, the assumptions are also made that the fringe shift 
is 0.02 too high at every zone and also that the fringe shift is 0.02 
too low at every zone. The density distributions calculated under these 
assumptions are not shown in figure 12 since they are almost indistin- 
guishable from each other and also from the density distribution 
obtained under the assumption of an uncertainty of 0.001. 

Equation (3) is now applied and is found to give, at zone 15, for 
AS = 0.02, a value of Ap of 0.00014 slug per cubic foot. The actual 
value of Ap is 0.00001 slug per cubic foot. For this case, therefore, 
equation (3) gives a very conservative estimate of the error in density. 

Applying equation (3) to a case where it will give the least con- 
servative estimate of error is desirable. The largest error in density 
occurs when the error in the fringe-shift curve is an oscillating one 
that changes sign at each zone. The density distribution obtained 
under this assumption of fringe-shift error is shown by the oscillating 
curve in figure 12 for |AS| = 0.02. At zone 15 the error in density 
is 0.0001 slug per cubic foot. This value is to be compared with the 
value 0.00014 slug per cubic foot calculated by equation (3). Equa- 
tion (3), therefore, appears to be a good approximation in this case, 
which 2.5 one that simulates those that occur in actual practice. 

Satisfying condition (e), that is, obtaining fringe shift with 
su i'i'i c i en i accuracy, is perhaps the most difficult requirement to meet 
in obtaining density distributions in three-dimensional flows. Because 
equation (3) so clearly illustrates the difficulty, the equation is 
used in another example to show how accurately the fringe shifts must 
be known for the density to be determined accurately. A typical posi- 
tion in figure 5(b) is chosen and the requirement is made that the 
density be obtained within 1 percent. A position is chosen about midway 
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along the contour on which the ratio of density to free-stream density 
is 2.0. For this position, 

_s 

X = 5.17 x 10 centimeters 
z = 2.3 X 10"^ centimeters 
k = O.II 67 cubic foot per slug 
i = 24 
A i;i = 0.11. 

p = 0.0070 slug per cubic foot 
Ap = 0.00007 slug per cubic foot 
Then, by equation (3), 


AS = 0.04 

In the present investigation, fringe shifts could generally be deter- 
mined with approximately twice this amount of uncertainty, or AS 
approximately equal to 0.1, for the interferogram of the flow at Mach 
number 1.62. For the interferogram taken at Mach number 1.30, howevex 
the fringe shifts were somewhat more uncertain than for the other inte. 
ferogram. The model is farther from the end of the nozzle, and turbu- 
lence in the mixing regions on the two sides of the jet through which 
the light passes causes appreciable irregularities in the fringe patter 

That the effect of shrinkage of the film during development was 
negligibly small was found by photographing a transparent scale. 
Furthermore, by enlarging photographs of a scale, it was found that dis- 
tortions introduced by the lens of the enlarger were also negligibly 
small . 

Condition (f) is that S(u) can be replaced by a straight line in 
each zone. Weyl has shown (reference 2) that this condition is suffi- 
ciently well satisfied everywhere except in zones in which the rate of 
change of S'(u) is large, that is, where the second derivative 
of S(u) has large absolute values. At shock waves the second deriva- 
tive of S(u) is very large. The result is that condition (f) is 
grossly violated and the calculated density overshoots at shock waves. 

At other places where the second derivative is large, but not so large 
as at shock waves, condition (f) can, by the use of smaller zone widths, 
generally be sufficiently well satisfied in that the error due to 
violation of condition (f) is acceptably small. 
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APPENDIX D 

COMPUTATIONAL PROCEDURE 


The negatives of the inter ferograms were about li inches in diame- 
ter. Prints were made in which the negatives were enlarged about 
16 diameters and these prints were used for all measurements. 

The computational procedure used is as follows: 

The axis of symmetry was drawn on each interferogram. On each no- 
flow interferogram the fringes were numbered arbitrarily. A cross sec- 
tion was drawn perpendicular to the axis of symmetry and the position 
of each fringe was measured along the cross section. Then fringe number 
was plotted as a function of position. A large scale that will permit 
reading the curve with the required accuracy had to be used. 

On the flow interferograms, lines perpendicular to the axis of sym- 
metry were drawn. These lines were spaced about l/lO inch apart and 
denoted the cross sections along each of which the fringe shift was to 
be measured. Again the fringes were arbitrarily numbered and the posi- 
tions of the fringes were measured along each cross section. The posi- 
tions of both the black and the white fringes were measured, because on 
part of the interferogram the fringes run almost parallel to the cross 
sections and the number of points is, therefore, less than is desirable. 

The next step was to determine the fringe shift at each measured 
point on each cross section. The position of each fringe is known. The 
position of the disturbed fringe was used in the previously mentioned 
plot of fringe number as a function of fringe position to determine the 
number of the undisturbed fringe at that position. Then the number of 
the undisturbed fringe was subtracted from the number of the disturbed 
fringe to give the fringe shift. 

Because the fringes had been arbitrarily numbered, all fringe 
shifts were adjusted by the amount necessary to make the undisturbed 
fringe at the same location in both flow and no-flow interferograms 
have the same number. This adjustment is tantamount to numbering the 
fringes on either the flow or the no-flow interferogram not arbitrarily 
but in such a way that the same number is given to the same position in 
the undisturbed portion of both interferograms. 

Because the free-stream density was to be used as the reference 
density in computing density distribution around the sphere and because 
the fringe shift must be assumed to be zero for the reference density. 
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the fringe shift in the free stream had to he determined and subtracted 
from all other fringe shifts. The fringe shift in the free stream was 
determined from the free-stream density which was calculated as described 
in appendix C. 

A cross section in the free stream was chosen as giving the posi- 
tions of the reference fringes. The fringes in the free stream were 
at a slight angle to the axis and, therefore, to the cross sections 
which were drawn perpendicular to the axis. This tilt of the free- 
stream fringes has an effect on the fringe shifts measured along each 
cross section in that it causes them to be different by an amount that 
is proportional to the tangent of the angle of tilt and to the distance 
of the cross section from the reference cross section. Proper adjust- 
ment of the measured fringe shifts along each cross section was made to 
eliminate the effect of the tilt. 

The next step was to plot fringe shift as a function of position. 

Here again, care must be taken to use a sufficiently large scale. The 
plot was then divided into zones, with the zeroeth zone at the axis of 
symmetry. Zone width is generally determined by a number of factors. 
Equation (3) shows that the greater the zone width the less the error 
in density due to violation of condition (e). On the other hand, the 
greater the zone width, the smaller the number of zones, the less well 
condition (f) is satisfied, and the smaller the number of experimental 
points that will be obtained. A compromise in zone width is, therefore, 
necessary. For the present investigation, the cross sections were 
separated into groups of nine each. The same zone width was used for 
all cross sections in a group and was such that the longest cross sec- 
tion in each group was divided into 50 zones. 

The next step was to read the fringe shift at each zone and then 
to determine, by subtracting, the change in fringe shift between adja- 
cent zones. This last quantity is the quantity - S^ + jl that is 

used in equation (2) for calculating v. The computations were made by 
the Bell computer. From the values of v the values of the ratio of 
the density to the free-stream density were computed. These values were 
used to obtain a plot of the density ratio against zone number for each 
cross section. A typical plot is shown in figure 13, which is for a 
cross section that is about three-fourths of the way between the nose 
of the shock wave and the nose of the sphere. 

The plot for each cross section has two regions where the inac- 
curacy or uncertainty is greater than for the rest of the curve. One 
of these regions is at and near the shock wave. There the inaccuracy 
is caused by failure to satisfy condition (f), and the result is over- 
shooting or undershooting at the shock wave. (in fig. 13 the shock wave 
is at zone ^3 and the correct density ratio behind the shock wave is 1.77*) 
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The density ratio across a shock wave is known, however, from well- 
established theory. The correct values of the density ratio behind the 
shock wave, therefore, can be shown on the plots of density ratio. The 
other region of relatively large inaccuracy is at and near the axis of 
symmetry (zones 0 to 5 in fig. 13). That the region of the axis is a 
region of relatively large inaccuracy can he seen from equation (3), 
which shows that the uncertainty in density is proportional to the coef- 
ficient A^ ^ . The coefficient ^ can he seen from table I to he 

larger for zone 0, which is at the axis of symmetry, than for any other 
zone. Unlike the other region of relatively large inaccuracy, however, 
for which the density ratio immediately behind the shock wave is known, 
the density-ratio distribution along the axis of symmetry is not known 
from theory. 
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TABLE I.- COEFFICIENTS A k * 



\l ( k + 1) 2 - l 2 - \/k 2 - 


2k + 1 


2 ] 


0 

1.00000 

1 

0.33333 

2 

.20000 

3 

.14286 

4 

.11111 

5 

.090909 

6 

.076923 

7 

. 066667 

8 

.058824 

9 

.052632 

10 

.047619 

11 

.043478 

12 

.040000 

13 

.037037 

14 

.034483 

15 

.032258 

16 

.030303 

17 

.028571 

18 

.027027 

19 

.025641 

20 

.024 390 

21 

.023256 

22 

.022222 

23 

.021277 

24 

.020408 

25 

.019608 

26 

.018868 

27 

.018182 

28 

.017544 

29 

.016949 

30 

.016393 

31 

.015873 

32 

.015385 

33 

.014925 

34 

.014493 

35 

.014085 

36 

.013699 

37 

.013333 

38 

.012987 

39 

.012658 

4o 

.012346 

4l 

.012048 

42 

.011765 

43 

.011494 

44 

.011236 

45 

.010989 

46 

.010753 

47 

.010526 

48 

.010309 

49 

.010101 

50 

.009901 

k / 


/ i 

0 


0.57735 

.21928 0.44721 
.14922 .17543 
.11400 .12427 
.092464 .097662 

.077856 .080873 
.067270 .069184 
.059236 .060529 
.052926 .053842 
.047837 .048509 

.043644 .044153 
.040129 .040523 
.037139 .037451 
.034565 .034816 
.032325 .032530 

.030359 .030528 
.028618 .028760 
.027067 .027186 
.025675 .025777 
.024419 .024507 

.023281 .023357 
.022244 .022311 
.021296 .021354 
.020425 .020477 

.019623 .019668 

.018882 .018922 
.018194 .018230 
.017555 .017587 
.016959 . 016988 
.016402 .016429 

.015881 .015905 
.015392 .015414 
.014932 .014952 
.014499 .014517 
.014090 .014107 

.013704 .013719 
.013338 .013352 
.012991 .013005 
.012662 .012675 
.012349 .012361 

.012052 .012062 
.011768 .011778 
.011497 .011506 
.011239 .011247 
.010992 .011000 

.010755 .010763 
.010529 .010536 
.010311 .010318 
.010103 .010109 
.009901 .009911 


1 2 


). 37796 

.15047 0.33333 
.10874 .13383 

.086800 .097879 
.072776 .078909 
.062887 .066709 
.055480 .058047 
.049696 .051514 

. 045041 .046381 
.041206 .042224 
.037988 .038781 
.035247 .035877 
.032880 .033390 

.030817 .031236 
.029001 .029349 
.027390 .027682 
.025950 .026198 
.024656 .024869 

.023486 .023669 
.022423 .022582 
.021452 .021592 

.020563 .020686 
.019745 .019854 

.018990 .019087 
.018291 .018377 
.017642 .017719 
.017037 .017107 
.016473 .016536 

.015945 .016003 
.015451 .015502 

.014986 .015033 
.014548 .014591 
.014135 .014175 

.013745 .013782 
.013376 .013410 
.013027 .013058 
.012695 .012724 
.012380 .012406 

.012080 .012105 
.011794 .011817 
.011522 .011543 
.011262 .011282 
.011013 .011032 

.010775 .010793 
.010547 .010564 
.010329 .010345 
.010120 .010134 
.009921 .009931 


3 4 


). 30151 

.12172 0.27735 
.089735 .11240 
.072842 .083335 
.061944 .067989 
.054176 .058074 

.048294 .050989 
.043652 .045610 
.039878 .041354 
.036739 .037883 
.034082 .034989 

.031800 .032532 
.029815 .030417 
.028072 .028573 
.026529 .026949 
.025150 . 02^y^Q 

.023912 .024219 
.022792 .023058 
.021776 .022006 
.020847 .021049 
.019996 .020175 

.019213 .019371 
.018490 .018631 
.017820 .017946 
.017198 .017311 
.016618 .016720 

.016077 .016169 
.015570 .015654 
.015094 .015171 
.014648 .014717 
.014226 .014290 

.013829 .013888 
.013453 .013507 
.013098 .013148 
.012761 .012807 
.012441 .012483 

.012137 .012176 
.011847 .011884 

.011571 .011605 

.011308 .011340 
.011056 .011086 

.010815 .010843 
.010585 .010611 
.010365 .010389 
.010153 .010176 
.009950 .009970 


5 6 


). 25820 

.10493 0.24254 
.078135 .098784 
.063993 .073802 

. 05^848 .060628 
.048306 .052107 
.043330 .046007 
.039384 .041360 
.036161 .037672 

.033468 .034654 
.031177 .032129 
.029200 .029978 
.027474 .028118 
.025951 .026492 

.024597 .025056 
.023383 .023777 
.022289 .022629 
.021296 .021592 
.020392 .020651 

.019563 .019792 
.018801 .019004 
.018098 .018279 
.017448 .017609 
.016843 .016989 

.016280 .016411 
.015754 .015873 
.015262 .015370 
.014801 .014899 
.014367 .014456 

.013958 .Ol4o4o 
.013572 .013648 
.013207 .013277 
.012862 .012926 
.012534 .012594 

.01222 3 .012279 
.011928 .011979 
.011646 .011694 
.011378 .011422 
.011121 .011163 

.010877 .010915 

.010643 .010679 
.010418 .010452 
.010204 .010236 
.0099901 .010030 


7 8 


). 22942 

.093603 0.21822 

,070117 .089160 
.057743 .066935 
.049740 .055235 
.044007 .047668 
.039637 .042247 

.036165 .038113 
.033321 .034826 
.030939 .032131 
.028907 .029872 
.027149 .027943 

.025609 .026273 

. 024248 .024809 
.023034 .023513 
.021943 .022356 
.020957 .021316 

.020061 .020375 
.019242 .019519 
.018490 .018736 
.017798 .018016 
.017158 .017353 

.016564 .016739 
.016011 .016169 
.015495 .015639 
.015013 .015143 
.014560 .014679 

.014135 .014244 
.013735 .013834 
.013357 .013449 
.013000 .013085 
.012662 .012740 

.012342 .012414 

.012038 .012105 
.011749 .011811 
.011473 .011531 
.011210 .011265 

.010960 .011010 
.010721 .010768 
.010492 .010536 
.010272 .010314 
.OIOO69 .010099 
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TABLE I.- COEFFICIENTS A k ^ - Continued 


11 

12 

13 

14 

15 

16 

17 

18 

0.20851 

.085295 0.20000 

.064150 .081893 0.19245 

.053027 .061686 .078868 0.18570 

.045836 .051065 . 0591+86 . 076155 0.17961 

.040683 .044199 .049305 .057506 .073704 0.17408 

.036752 .039281 .042727 .047716 .0557H .071476 0.16903 

.033625 .035527 .038014 .041392 .046271 .054074 .069438 0.16440 


19 

.031061 

.032541 

.034417 

.036862 

.040176 

.044950 

.052574 

.067565 0.16013 

20 

.028909 

.030090 

.031555 

.033405 

.035809 

.039060 

.043736 

.051192 

.065036 

21 

.027070 

.028033 

.029205 

.030653 

.032478 

.03^2 

.038033 

.042615 

.049913 

22 

.025477 

.026274 

.027232 

.028394 

.029825 

.031623 

.033949 

.037082 

.041576 

23 

.024079 

.024749 

.025544 

.026496 

.027647 

.029060 

.030833 

.033121 

.036200 

24 

.022841 

.023410 

.024080 

.024872 

.025817 

.026956 

.028352 

.030099 

.032351 

25 

.021735 

.022224 

.022794 

.023462 

.024250 

.025187 

.026314 

.027692 

.029415 

26 

.020740 

.021163 

.021654 

.022224 

.022890 

.023673 

•024601 

.025716 

.027077 

27 

.019839 

.020208 

.020634 

.021126 

.021695 

.022358 

.023135 

.024055 

.025157 

28 

.019018 

0 019343 

.019715 

.020143 

.020635 

.021202 

.021861 

.022632 

.023543 

29 

.018267 

.018554 

.018882 

.019257 

.019685 

.020176 

.020741 

.021396 

.022161 

30 

.017577 

.017832 

.018123 

.018453 

.018829 

.019257 

.019747 

.020309 

.020959 

31 

.016940 

.017168 

.017427 

.017720 

.018052 

.018428 

.018856 

.019344 

.019903 

32 

.016350 

.016555 

.OI6786 

.017048 

.017343 

.017676 

.018052 

.018479 

.018964 

33 

.015802 

.015987 

.016195 

.016429 

.016693 

.016989 

.017322 

.017698 

.018123 

34 

.015291 

.015458 

.015646 

.015857 

.016094 

.016359 

.016656 

.016989 

.017364 

35 

.014814 

.014966 

.015136 

.015327 

.015540 

.015778 

.016044 

.016342 

.016675 

36 

.014367 

.014505 

.014660 

.014833 

.015026 

.015241 

.015481 

.015747 

.016045 

37 

.013947 

.014073 

.014215 

.014373 

.014548 

.014743 

.014959 

.015199 

.015466 

38 

.013552 

.013668 

.013798 

.013942 

.014102 

.014279 

.014475 

.014692 

.014932 

39 

.013180 

.013286 

.013405 

.013537 

.013683 

.013845 

.014024 

.014221 

.014439 

1+0 

.012828 

.012926 

.013036 

.013157 

.013291 

.013439 

.013602 

.013782 

.013980 

1+1 

.012495 

.012586 

.012687 

.012799 

.012922 

.013058 

.013207 

.013372 

.013552 

1+2 

.012180 

.012264 

.012357 

.012460 

.018574 

.012699 

.012837 

.012987 

.013152 

43 

.011880 

.011958 

.012045 

.012140 

.012246 

.012361 

.012487 

.012626 

.012778 

44 

.011596 

.C11668 

.011749 

.011837 

.011934 

.012041 

.012158 

.012286 

.012426 

45 

.011325 

.011392 

.011467 

.011549 

.011640 

.011739 

.011847 

.011965 

.012094 

46 

.011067 

.011130 

.011199 

.011276 

.011360 

.011452 

.011552 

.011662 

.011781 

47 

.010821 

.010879 

.010944 

.011016 

.011094 

.011180 

.011273 

.011375 

.011485 

48 

.010585 

.oio64o 

.010701 

.010768 

.010841 

.010921 

.011008 

.011102 

.011205 

49 

.010360 

.010412 

.010468 

.010531 

.010599 

.010674 

.010755 

.010843 

.010939 

50 

.010149 

.010188 

.010248 

.010307 

.010366 

.010436 

.010515 

.010594 

.010693 

yi 

11 

12 

13 

l4 

15 

16 

17 

18 

19 


.028776 

.026501 

.024633 

.023062 

.021717 

. 02051+8 

.019520 

.018607 

.017788 

.01701+9 

.016377 

.015763 

.015199 

.011+679 

.011+198 

.013751 

.013334 

.01291+1+ 

.012578 

.012235 

.0119U 

.011605 

.011316 

.01101+3 

.010782 

20 
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TABLE I.- COEFFICIENTS A kj ± - Continued 


21 

0.15250 










22 

.062742 0.14907 









23 

.047620 

.061350 0.14586 








24 

.039708 

.046585 

.060047 0.14286 







25 

.034609 

.038864 

.045616 

.058824 0.14003 






26 

.030959 

.033888 

.038071 

.044704 

.057672 0.13736 





27 

.028176 

.030328 

.033211 

.037325 

.043845 

.056585 0.13484 




28 

.025960 

.027613 

.029733 

.032572 

.036621 

.043034 

.055558 0.13245 



29 

.024141 

.025452 

.027082 

.029173 

.031969 

.035956 

.042266 

.054585 0.13019 


30 

.022610 

.023676 

.024971 

.026581 

.028642 

.031399 

.035325 

.041538 

.053661 0.12804 

31 

.021300 

.022184 

.023238 

.024517 

.026106 

.028140 

.030857 

.034727 

.0^846 

.052782 

32 

.020161 

.020905 

.021780 

.022823 

.024087 

.025656 

.027663 

.030343 

.034158 

.040188 

33 

.019159 

.019794 

.020532 

.021398 

.022430 

.023679 

.025229 

.027210 

.029854 

.033616 

34 

.018268 

.018817 

.019447 

.020178 

.021036 

.022056 

.023291 

.024822 

.026778 

.029388 

35 

.017470 

.017948 

.018492 

.019117 

.019842 

.020691 

.021700 

.022921 

.024434 

.026366 

36 

.016749 

.017169 

.017644 

.018184 

.018804 

.019522 

.020362 

.021361 

.022569 

. 024064 

37 

.016095 

. 016466 

.016883 

.017354 

.017891 

.018506 

.019217 

.020049 

.021037 

.022232 

38 

.015496 

.015826 

.016196 

.016611 

.017079 

.017611 

.018221 

.018926 

.019750 

.020728 

39 

.014946 

.015242 

.015571 

.015939 

.016351 

.016816 

.017345 

.017949 

.018647 

.019464 

4o 

.014439 

.014705 

.015000 

.015328 

.015693 

.016103 

.016565 

.017090 

.017689 

.018381 

4l 

.013969 

.014210 

.014475 

.014769 

.015095 

.015459 

.015867 

.016325 

.016845 

.017440 

42 

.013532 

.013751 

.013991 

.014256 

.014549 

.014873 

.015235 

.015640 

.016095 

.016611 

43 

.013125 

.013324 

.013543 

.013782 

.014046 

.014338 

.014661 

.015021 

.015423 

.015875 

44 

.012745 

.012926 

.013125 

.013343 

.013582 

.013846 

.014136 

.014457 

.014815 

.015214 

45 

.012387 

.012554 

.012736 

.012935 

.013153 

.013391 

.013653 

.013942 

.014262 

.014617 

46 

.012052 

.012205 

.012372 

.012554 

.012753 

.012970 

.013208 

.013469 

.013756 

.014074 

47 

.011736 

.011877 

.012031 

.012198 

.012380 

.012578 

.012795 

.013032 

.013292 

.013578 

k 8 

.011437 

.011568 

.011710 

.011864 

.012031 

.012213 

.012411 

.012626 

.012862 

.013121 

49 

.011155 

.011276 

.011407 

.011549 

.011704 

.011871 

.012052 

.012249 

.012464 

.012700 

50 

.010891 

.011000 

.011119 

.011257 

.011396 

.011554 

.011723 

.011891 

.012089 

.012307 

k / 
/i 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 
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TABLE I.- COEFFICIENTS i - Continued 


31 

0.12599 










32 

.051945 

0.12404 









33 

.039560 

.051147 0.12217 








3^ 

.033099 

.038961 

.050385 0.12039 







35 

.028943 

.032605 

.038389 

.049655 0.11868 






36 

.025973 

.028517 

.032133 

.037841 

.0148957 0.11704 





37 

.023710 

.025596 

.028110 

.031681 

.037315 

.048287 0.11547 




38 

.021910 

.023371 

.025236 

.027720 

.031247 

.036811 

.047644 0.11398 



39 

.020432 

.021601 

.023047 

.024890 

.027345 

.030830 

.036327 

.047026 0.11251 


4o 

.019190 

.020148 

.021305 

.022735 

.024559 

.026985 

.030430 

.035862 

.046431 0. 11111 

41 

.018126 

.018927 

.019876 

.021021 

.022436 

.024240 

.026640 

.030045 

.035414 

.045858 

42 

.017201 

.017881 

.018675 

.019614 

.020748 

.022149 

.023933 

.026307 

.029674 

.034982 

43 

.016387 

.016972 

.017646 

.018432 

.019363 

.020486 

.021872 

.023637 

.025986 

.029317 

44 

.015663 

.016171 

.016751 

.017419 

.018199 

.019121 

.020233 

.021605 

.023353 

.025677 

45 

.015014 

.015460 

.015964 

.016539 

.017202 

.017974 

.018888 

.019989 

.021348 

.023078 

46 

.014427 

.014822 

.015264 

.015764 

.016335 

.016992 

.017758 

.018663 

.019754 

.021100 

47 

.013894 

.014245 

.014636 

.015076 

.015572 

.016138 

.016790 

.017549 

.018446 

.019527 

48 

.013406 

.013720 

.014069 

.014458 

.014894 

.015387 

.015948 

.016594 

.017347 

.018236 

49 

.012957 

.013240 

.013553 

.013899 

.014286 

.014719 

.015208 

.015765 

.016406 

.017152 

50 

.012545 

.012792 

.013079 

.013386 

.013733 

.014119 

.014554 

.015035 

.015584 

. 016224 


31 


32 33 


34 35 3 6 37 38 


39 


4o 
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TABLE I.- COEFFICIENTS A kjl - Concluded 


4 l 

0.10976 









42 

. 045306 0.10847 








^3 

.034566 

.044774 0.10721 







44 

.028972 

.034164 

.044260 0.10600 






45 

.025378 

.028639 

.033776 

.043763 0.10483 





46 

.022813 

.025090 

.028318 

.033401 

.043282 0.10370 




bl 

.020860 

.022557 

.021(811 

. 028007 

.033038 

.042817 

0.10260 



48 
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(a) Full field of view. 
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Figure 1.- Undisturbed fringes for analyzing interferogram of flow at 

Mach number 1.30. 
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Figure 1.- Concluded 
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(a) Full field of view. 



L-6J4.897 


Figure 2.- Interferogram of flow around sphere at Mach number 1-30. 
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(b) Region used in analysis . 
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Figure 2.- Concluded 
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(a) Full field of view. 

Figure 3-- Undisturbed fringes for analyzing interferogram of flow at 

Mach number 1.62. 
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(b) Region used in analysis. 
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Figure 3.- Concluded 
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(a) Full field, of view. 


L-61|901 

Figure 4.- Interferogram of flow around sphere at Mach number 1.62. 
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Figure 4.- Concluded. 

















(a) M = 1.30. 

Figure 5 .- Contours of constant ratio of density to free-stream density. 
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Figure 5 «- Concluded. 
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Figure 6.- Variation along axis of ratio p/pq* 


46 


NACA TN 2173 



Shock Sphere 

Percent of distance from shock wave to sphere 


Figure 7>- Variation along axis of ratio p/p s . 
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Figure 8.- Variation along axis of local Mach number. 
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Figure 9 .- Variation of the density ratio p/pq with angular distance 
from the stagnation point along the surface of the sphere. M = 1.62. 
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Figure 10.- Concentric zones in disturbance. 
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Figure 11.- Refraction of light ray by sphere. 
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Figure 12.- Assumed distribution of density, corresponding distribution 
of fringe shift, and calculated distributions of density. 


NACA-Langley - 8-23-50 -925 



VJ1 

ro 


Figure 13*- Variation of ratio 


P^Pq along typical cross section. 
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